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NUMERICAL SOLUTION OF THE STEADY- STATE NAVIER- STOKES EQUATIONS 
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and 

James N. Moss 
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ABSTRACT 

The steady- state Navier- Stokes equations are solved for 
hypersonic flow about blunt ax i symmetric bodies, The equations 
of motion are solved by successive approximations using an 
implicit finite-difference scheme. The results are compared 
with viscous shock- layer theory, experimental data, and time- 
dependent solutions of the Navier-Stokes equations. It is demon- 
strated that viscous shock- layer theory is sufficiently accurate 
for the range of flight conditions normally encountered by entry 
vehicles. 

INTRODUCTION 

The solution for the flow field about blunt ax i symmetric 
bodies at zero angle of attack has been determined for moderate 
to high Reynolds number conditions using the combination of 
inviscid and boundary- layer theories or viscous shock- layer 
theory. However, in a recent publication, Li [1] questions the 
validity of existing theories in the range of Reynolds numbers 
encountered by entry vehicles in the altitude range 200,000 to 
300,000 feet. Li presents a time- dependent finite-difference 
solution of the Navier-Stokes equations and applies the method 
to typical shuttle trajectory points in this altitude range. 

This method of solution is severely restricted by storage 
requirements and excessive computing cost and, as a result, 
limited in application to the nose region of the body. 

In a previous publication, Davis [2] presents a solution of 
the viscous shock- layer equations and shows good agreement -with 
experimental data for Reynolds numbers as low as 22. The viscous 
shock-layer equations presented by Davis [2] contain all terms up 
to 0 (e*), where e is the Van Dyke [33 expansion parameter 1 
defined as 

€ = M T ref )/p„ u „ V l/2 

and 


L-9200 



T . = U 2 /C 
ref 00 P 

= body nose radius 


U M ~ free- stream velocity 


[x = coefficient of viscosity 


Consideration of the order of magnitude analysis of the 
Navi er- Stokes equations given by Van Dyke [3l and the viscous 
shock- layer equations presented by Davis [2], shows that the 
solution of the two sets of equations will have significant dif- 
ferences only when e ~ 0(l). This suggests that a method of 
successive approximations in which the starting solution is given 
by the viscous shock- layer equations can be used to determine the 
solution of the steady-state Navier-Stokes equations for most 
flight conditions of practical interest. 


The solutions of the Navier-Stokes equations presented herein 
are restricted to a perfect gas having a constant ratio of spe- 
cific heats and the bow shock is assumed distinct. However, 
appropriate slip and jump boundary conditions are included. The 
solutions of the Navier-Stokes equations are compared with viscous 
shock- layer theory, experimental data, and the transient solution 
presented by Li Cl]. 


METHOD OF SOLUTION 

The Navier-Stokes equations are expressed in the shock-body 
oriented coordinate system defined by Davis [2]. After deter- 
mining the solution of the viscous shock- layer equations using 
the implicit finite-difference scheme described in Reference [2], 
the terms of higher order are evaluated using these flow-field 
data. The higher order terms are held constant during the solu- 
tion for the first approximation to the Navier-Stokes equation. 
After the flow- field data for the first approximation are deter- 
mined, the higher order terms are reevaluated and held constant 
during the solution for the second approximation. This procedure 
is repeated until the flow-field data corresponding to successive 
approximations converge within a specified limit. For details of 
the numerical method, the reader is referred to the paper pre- 
sented by Davis [2], 

! RESULTS AND DISCUSSION 

i 

A comparison of the temperature profiles along the stagna- 
tion streamline and at a distance of 0.7 nose radii downstream 
of the Stagnation point is shown in Figure 1 for solutions corre- 
sponding to the steady- state Navier-Stokes equations, the viscous 
shock- layer equations, and Li's transient solution of the Navier- 
Stokes equations. The flight conditions for this case correspond 
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Figure 1. Comparisons of temperature profiles 
for transient and steady- state solutions. 

to a Mach number of 2h at an altitude of 250,000 feet, a free- 
st ream Reynolds number of 1^,^00 based on a 2-foot nose radius, 
and an e of O.Ot. The body considered is a 20° hyperboloid. 

For these conditions, the present solution of the Navier-Stokes 
equations and the solution of the viscous shock-layer equations 
differed by less than 0.5$ and these results are shown as a 
single curve. The solutions given by Li are seen to exhibit 
large oscillations at both body locations. Since the maximum 
number of node points across the shock layer in Li's solution is 
less than 35> these oscillations can probably be removed by 
increasing the number of node points. For the present analysis, 
a minimum of 50 grid points was necessary for adequate resolu- 
tion of this problem. As indicated by Li, this is apparently a 
prohibitive requirement in his solution. It is noted that the 
temperature behind the shock at the downstream station computed 
using the transient method of solution is approximately 15$ lower 
than predicted by the present method. The temperature predicted 
by Li closely corresponds with the value obtained assuming that 
the local shock angle is the same as that at the body. The com- 
parison of the present solution of the steady- state Navier-Stokes 
equations and the viscous shock-layer equations with the tran- 
sient solution of Reference [1] indicates that the transient 
solution is an unacceptable alternate method for solving the 
viscous blunt-body problem. 

The present solution of the Navier-Stokes equations is com- 
pared with the solution of the viscous shock- layer equations i and 
experimental data for Reynolds numbers in the range 90 to 31, 160 
in Figure 2. The experimental data are given by Little [4]. The 
comparisons of the total drag coefficient, C^, shown in Figure 2 
are representative of the differences noted in heat transfer, 
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Figure 2. Total drag coefficient for a 1(5° hyperboloid. 

skin friction, and surface pressure distribution using the two 
sets of governing equations. For the higher Reynolds number 
case, the two solutions differed by less than 0.1# and are shown 
as a single curve. At a Reynolds number of 1035, the maximum 
difference m the two solutions is less than 5#, and at a 
Reynolds number of 90, the two solutions differ by approximately 
15 i V n the stagnation region and less than 2 % in the downstream 

anJPtu' w Th ? so ^ t 1 ion of both the viscous shock-layer equations 
and the Navier-Stokes equations is in good agreement with the 

experimental data. Figure 3 shows the skin friction and heat- 
r nsfer distributions for the case of a Reynolds number of 90. 

The solutions of the Navier-Stokes equations presented 
herein are seen to reduce to the solution of the viscous shock- 
yer equations for high Reynolds number flow and significant 
erences occur only at very low Reynolds numbers. These 
r6 f U 3 are ^- n agreement with the requirements of Van Dvke's f3l 
order °f magnitude analysis. The cases presented were solved to 
a distance of 8 nose radii downstream of the stagnation point 
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Figure 3. Comparisons of results for a 45° hyperboloid. 

using a CDC-6600 digital computer. The maximum computing time 
was less than 3 minutes and the storage requirement is less 
than 50 ,,OOOq. 

The results of this investigation reveal that there is little 
need for a solution of the Navier- Stokes equations for viscous 
blunt-body flows in the range of flow conditions of practical 
interest. However, a method which is accurate and computationally 
efficient has been developed for their solution. The numerical 
method with proper modifications can be extended to merged layer 
problems and reacting gas mixtures. 
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